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The onset of undamped wave propagation in noisy self-oscillatory media is identified with 
a Hopf bifurcation of the corresponding effective dynamical system obtained by properly 
renormalizing the effects of noise. We illustrate this fact on a dense array of nonlocally 
coupled phase oscillators for which a mean-field idea works exactly in deriving such effective 
dynamical equations. 

Unlike conservative oscillatory media, dissipative self-oscillatory media admit 
traveling waves without decay even in the presence of noise or other sources of ran- 
domness. This ability of wave transmission in random self-oscillatory media should 
be functionally relevant in a variety of living organisms for which randomness is 
unavoidable. For this type of systems, a critical strength of randomness is gener- 
ally expected to exist such that below which the system is capable of sustaining 
undamped traveling waves. This critical point should be identical with the Hopf 
bifurcation point of an effective dynamical system obtained by properly renormaliz- 
ing the effects of noise. It is reasonable to expect that a theory could be developed 
unambiguously on this issue in the particular case when each oscillator couples with 
sufficiently many oscillators, because a mean-field idea should be applicable then. 
In the present article, we carry out this program for nonlocally coupled phase os- 
cillators with noise, and show how an effective dynamical equation can be derived, 
and how it is reduced to a small-amplitude equation near the bifurcation point. Our 
theory may be regarded as a natural extension of a previous theory El) on the onset 
of collective oscillation for globally coupled phase oscillators with noise. 

Imagine an infinitely long array of nonlocally coupled oscillators which are dis- 
tributed densely and subject to additive noise. By taking the continuum limit, the 
phases (j){x, t) of the oscillators are assumed to obey the following Langevin-type 
equation. 

(/.(x, t)=uj+ G{x- x')r(0(x, t) - (l){x', t')) + t). (1) 



Here oj is the natural frequency, the second term represents the nonlocal coupling, 
and the last term gives additive noise. The phase coupling function T depends 
only on the phase difference and satisfies the in-phase condition r'{0) < OEl-* , while 
its strength G depends on the distance. Specifically, we assume the form G{x) = 
7exp(— 7|x|)/2 whose integral is normalized. The noise is assumed to be white 
Gaussian with vanishing mean, i.e., < (,{x,t) >= and < (,{x,t)^{x' ,t') >= 2D5(t — 
t')5{x — x'). Since we are working with an oscillator continuum, the coupling radius 
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7^^ is so large that infinitely many local oscillators may fall within it. This allows us 
to apply a mean field idea similar to the one applied successfully to globally coupled 
oscillators involving randomness. We also remark that the model given by Eq. (1) is 
not merely our invention but the one derived under suitable conditions from a certain 
class of reaction-diffusion systems after eliminating adiabatically a highly diffusive 
chemical componentH^. 

Periodic traveling waves in our oscillatory field correspond to a family of solutions 
of Eq. (1) having the form = kx+Qt possibly disturbed by the noise. In the absence 
of noise, it is easy to confirm their existence and analyze their stability for general 
phase-coupling function. In the simple case given by 

r(</>) = - sin((/) + a) (|a|<7r/2), (2) 

for instance, the traveling waves of sufficiently small k turn out always stable. The 
additive noise is now switched on. Numerical simulation for \a\ not too close to 7r/2 
shows that as the noise becomes stronger, traveling waves of shorter wavelengths 
(still satisfying the stability condition for the noiseless case) successively become 
unstable. For sufficiently strong noise, no traveling waves can exist, implying that 
the medium is no longer self-oscillatory because the traveling wave with infinite 
wavelength corresponds to a uniform oscillation. In order to explain these results, 
we will now convert Eq. (1) to a nonlinear Fokker-Planck equation and apply the 
center-manifold reduction to it. 

In terms of the number density n('0, x, t) = 5{(j){x, t) Eq. (1) can be rewrit- 
ten as 

d(p{x,t) r dx'G{x-x') di;'r{<p{x,t) - i;')n{'^' ,x' ,t) + C{x,t). 



dt J~oo Jo 

= V{(l)ix,t),x,t)+Cix,t). (3) 

The drift velocity V introduced above is expressed as 

V{ip,x,t)=uj+ dip'Z{'>P',x,t)r{tl; -ip')di/j', (4) 







where 

Z{il}',x,t)= I dx'G{x-x')n{iP\x',t). (5) 



The quantity Z may be interpreted as a space-time dependent mean field. If we 
regard Z as an externally given quantity, Eq. (3) represents a non-autonomous one- 
oscillator equation. Thus, the normalized probability density I{ip,x^t) such that 
(j){x,t) takes value ip is governed by the Fokker-Planck equation cP 

'Ji^^-^[vi,,..m,,..t)],D'-^^ (0) 

We now take a statistical average of the dynamical variable n{ip',x',t) appearing 
in the mean field Z. This is allowed because Z is given by a weighted sum over 
infinitely many oscillators which are individually subject to random forcing so that 
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the statistical fluctuation of Z may safely be neglected. Noting that the average 
< n{ilj',x',t) > is identical with f{ij;',x',t) by definition, we finally find a nonlinear 
evolution equation for / in the form 



df{ij,x,t) ^ d_ 

dt dip 



oo i'2n 



dx'G{x-x') / d'ip'r{i;-ij')f{ij',x',t)}f{'ip,x,t) 

oo Jo 



Wave propagation and collective dynamics could be discussed on the basis of this 
deterministic kinetic equation. In particular, as is shown below, one may apply 
the center-manifold reduction near the critical noise strength D = Dc at which the 
trivial constant solution /o = (27r)~^ loses stability. For this purpose, it is conve- 
nient to rewrite Eq. (7) in terms of the fluctuation p{'ip,x,t) defined by f{ijj,x,t) = 
fo + pii'lXit). The 27r-periodic functions p{ip) and -r(V') are further developed into 
Fourier series like p{il}, x, t) = 27r~-'^ Z^z^o Pii^-: exp(i^^/;) and r{ip) = eyip{ilip). 
Using these expressions, Eq. (7) becomes 



/oo 
dx'G{x - x')pi{x',t) 
-oo 

POO 

-il / dx'G{x-x') ^ rmPmix,t)pi_rn{x,t). (8) 



dpi{x,t) 
dt 



Regarding the loss of stability of the constant solution /o, spatially uniform per- 
turbation is the most dangerous. This is because by definition the uniform Fourier 
component of G is the largest of all its Fourier components, thus making the linear 
growth rate of the uniform perturbation of pi the largest. Therefore, the sign of this 
growth rate ai defined by 

= -fD + l%ri) (9) 

determines the stability of /o with respect to the fluctuation pi. Note that '^{F^i) = 
—Q{ri), so that (7/ is independent of the sign of /. Suppose that ai as a function of 
I takes the largest value at I = it A (A > 0). As the noise becomes weaker, a Hopf 
bifurcation occurs at D = Dc = '<s{rx)/X, the corresponding frequency being given 

hy n=-i\{u; + ^{rx) + ro). 

The center-manifold reduction can be achieved almost in parallel with the case 
of globally coupled phase oscillators with noise . The only difference is the existence 
of spatial integrals in the governing equation. Since the characteristic spatial scale is 
expected to be very large near the bifurcation point, which can actually be confirmed 
from the results obtained, one may develop p(x') into a Taylor series like p(x') = 
p{x) + {x' — x)p'{x) + {x' — x)^p"{x)/2 + - ■ ■ in Eq. (8), and neglect higher order effects. 
It turns out that the replacement p{x') with p{x) is allowed in the nonlinear terms, 
while the second order derivative of p{x) must be included in the linear term. In 
this way, Eq. (8) takes the form of a reaction-diffusion system, and its reduced form 
near Dc is the complex Ginzburg-Laudau equation. As usual, a complex amplitude 
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A{x^t) related to p{'ip,x,t) is introduced as 

p{i), x,t) = — (a{x, t) exp[i(AV' + f2t)] + c.c.) . (10) 

By applying the standard method, the reduced equation near D = becomes 

f)A ff' A 

— = X\D^-D)A + d^-g\A\^A, (11) 
where d and g are complex coefficients given by 

7^ ' ^ 29(A)-i5?(A)+zr2A- ^ ^ 

We now consider the simple coupling function given by Eq. (2), and show some 
results for this particular case. Obviously, we have A = 1, I±i = ±iexp(ia)/2 and 
all /] with 1^1 > 2 vanish. The quantities d and g are simplified as 

d = = ^ exp(iQ!), (13) 

' ' (14) 



2Si(ri) -i5R(ri) 4 cos a + 2i sin a 

In this particular case, the bifurcation is supercritical because ^{g) > 0. 

From the well-known properties of the complex Ginzburg-Landau equation, Eq. 
(11) admits a family of traveling plane waves Ak{x, t) in the form A^ = Rk exp[z(fcx + 
i7fci)]. The stability condition of the uniform solution Aq is given by 1 + ciC2 > 0, 
where ci = S(d)/K(ci) and C2 = Q{g)/^{g). For the sine-coupling given by Eq. 
(2), the corresponding condition becomes tan'^ a < 2 or \a\ < 0.9553.... Under this 
condition, plane waves with k up to some critical value of order (D^ — D)^/'^ can 
propagate through the medium stably. This is consistent with the aforementioned 
results of our numerical simulation on the stochastic system Eq. (1). 

Note that if |a| exceeds the above critical value, the uniform oscillation becomes 
unstable and phase turbulence sets in. In terms of the original model given Eq. (1), 
this means that the noise can strongly be amplified to a turbulent level. Note also 
that Eq. (11) admits the Nozaki-Bekki hole solutions & for which the existence of 
amplitude degree of freedom is crucial. If the phase model Eq. (1) which apparently 
lacks the amplitude variables can also exhibit hole structures, this should solely be 
due to a strong phase scattering caused by the noise in a local group of phase os- 
cillators by which the effective amplitude there becomes vanishingly small. Detailed 
numerical study associated with these remarks will be reported elsewhere. 
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